The path-coalescence transition and its applications 
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We analyse the motion of a system of particles subjected a random force fluctuating in both space 
and time, and experiencing viscous damping. When the damping exceeds a certain threshold, the 
system undergoes a phase transition: the particle trajectories coalesce. We analyse this transition 
by mapping it to a Kramers problem which we solve exactly. In the limit of weak random force we 
characterise the dynamics by computing the rate at which caustics are crossed, and the statistics of 
the particle density in the coalescing phase. Last but not least we describe possible realisations of 
the effect, ranging from trajectories of raindrops on glass surfaces to animal migration patterns. 



This paper discusses a surprising phase transition il- 
lustrated in Fig. n This figure shows trajectories for a 
system of particles with positions Xi(t) subjected to a 
random force fluctuating in both space and time, and 
experiencing a resistive force proportional to their ve- 
locities (the equations of motion are defined by 0) and 
(J2J) below). The motion of any one particle is obviously 
diffusive. Two particles with very close positions and 
momenta must follow similar trajectories, at least for a 
while. Diffusive motion usually tends to reduce inhomo- 
gencitics in density, and we might expect that the motion 
should resemble the simulation in Fig. However, a 
slight increase in the resistive force causes a phase tran- 
sition to the situation shown in Fig. which we term 
'path coalescence'. This effect has been described in a 
paper by Deutschi, who gave a theoretical analysis and 
numerical simulations indicating the existence of a phase 
transition in this model. 

Here we point out that these and related equations of 
motion have a very broad range of applications in the 
physical sciences, ranging from tracks of raindrops on 
randomly contaminated glass surfaces, to energetic elec- 
trons in disordered solids. In the 'overdamped' limit, 
where inertia is negligible, the model could describe the 
response of animals to random fluctuations of their envi- 
ronment. Possible applications of this type include mi- 
gration tracks of mammals and clustering of microorgan- 
isms. Because of the large variety of applications we 
believe that the path-coalescence effect deserves to be 
thoroughly understood. This paper makes the follow- 
ing contributions. First, adopting an approach^ used in 
the theory of Anderson localisation, we map the equa- 
tions of motion to a Kramers problem which we solve 
exactly. This enables us to obtain an exact criterion for 
the phase transition. Second, wc characterise the dy- 
namics by computing the rate at which trajectories (in 
Fig. ^) cross caustics. Third, in the coalescing phase, 
we determine the statistics of the particle density by cal- 
culating its pair-correlation function. This allows us to 
deduce, at time t, the expected number of particles con- 
densing into a trail. Fourth we argue that the model we 
consider (which is more general than that put forward 
ini) exhibits a complex phase diagram. Finally we con- 
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FIG. 1: Positions of 25 particles with viscous dissipation 
subjected to a spatially and temporally varying random force. 
The dynamics are described by eqs. 0,0. In both cases 
m = t = 1, £ = 0.05, e = 0.001; aj^ 0.04, b 7 « 0.1. 



elude with a discussion of some of the applications of the 
path-coalescence effect. Given the ubiquity of diffusive 
dynamics, the effect we describe here is bound to occur 
in a wide variety of different contexts. 

We consider trajectories for a system of independent 
particles with positions Xi(t), and momenta Pi(t). The 
equations of motion for any particle are 

x=p/m, p = f(x,t) -7p (1) 

where x = dx/dt, and 7 characterises the strength of 
the viscous damping. The statistical properties of the 
force f(x,t) are translationally invariant in both space 
and time, 

(f(x,t))=0, (f(x,t)f(x',t'))=c(x~x',t-t'). (2) 

(A) is the ensemble average of A. The correlation func- 
tion decays rapidly as \x — x'\ — > 00 and as \t — t'\ — > 00. 
A suitable form (adopted in Fig. ^) is c(Ax,At) = 
e 2 exp(— Ax 2 /2{; 2 ) exp(— Ai 2 /2r 2 ) where e denotes the 
magnitude of the force. The dynamics of the model 
is characterised by two independent dimensionless vari- 
ables: x = e7 " 2 / TO £ is a dimensionless measure of the 
strength of the random force, and the motion is over- 
damped when i/E7T> 1 . 
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Phase transition. The phase transition is determined 
by the fraction of initial conditions for which the separa- 
tion a pair of infinitesimally close trajectories approaches 
zero as t — ► oo: this is either or 1. The separation Sx 
of two nearby trajectories has a lognormal distribution. 
We have 



(]n\6x(t)/6x(0)\) = Xt, 



(3) 



we call A the Liapunov exponent. When A < 0, almost 
all nearby trajectories eventually merge, conversely when 
A > nearby trajectories almost certainly diverge. The 
condition for the phase transition is therefore A = 0, and 
evaluation of A is also valuable because it gives the rate 
of coalescence. Linearising the equations of motion JJJ 
gives Sx = Sp/m, and Sp = —jSp + d x f(x,t) Sx. Here Sx 
and Sp are small separations between pairs of trajectories. 
As t — > oo the ratio between Sp and Sx has a stationary 
distribution. We therefore write X = Sp/Sx, and find 
equations of motion in terms of X and Sx: 

Sx = XSx/m (4) 
X = d x f(x(t),t)- 1 X-X 2 /m. (5) 

Since the distribution of X remains stationary, eq. Q 
implies that 



A = (X)/m. 



(6) 



Consider the dynamics of X. When —X is sufficiently 
large, the noise term 9 x /(x(t),t) may be neglected and 
(f5j) implies that X reaches — oo in finite time. This point 
is a caustic, where Sx passes through zero and X jumps 
from — oo to +oo. As t — > oo, one obtains a stationary 
distribution, with X going to — oo at a rate — J. 

We now specialise to the case where the displacement 
of X during the correlation time r is small compared to 
the deterministic terms (x <C 1), adopting the approach 
used in£. In this limit the probability density P(X, t) for 
X satisfies a Fokker-Planck equation 3 - 



8 t P = d x [(X 2 /m + 1 X)P] 



Vd x P. 



with diffusion constant 
1 d 2 



V 



2dx 2 



dt c(x, t) 



(7) 



(8) 



We require a steady-state solution of eq. (JHJ), P(X), with 
a constant flux J, satisfying J = v(X)P — VdP/dX with 
v(X) = —X 2 /m — 7A. It is convenient to express the 
solution in terms of dimcnsionless variables: 



P(X) = a((Z,n) , z = 



X 



1V) 1 / 3 



si = 



7?n 



2/3 



OV3 



(9) 



where the constant a and function £ are to be deter- 
mined. The exact stationary solution of JJJ is 



C(Z,fl) =ex-p[-(f>(Z,n)] I dY cxp[0(Y,ST)] 

J — OO 



(10) 




FIG. 2: a Liapunov exponent: theory (line) and numeri- 
cal experiments for x = 0.1 (o), 0.02 (□), and 0.001 (o). b 
Rate |J| of crossing of caustics, theory (line), and numerical 
experiments for \ = 0.02 (□). 



Here (j>(Z, fi) is proportional to the integral of v(X) 
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The flux J is determined by normalisation. Using @ we 
obtain the Liapunov exponent 



A = (P/m 2 ) 1/,3 *i(0)/*o(0) 



(12) 



where #„(fi) = ^dZ Z n ({Z,Q,). The Liapunov ex- 
ponent is shown in Fig. [2Jt. Consider the form of HUH . 
When f2 is large, the integral is approximately constant 
over the interval between the maximum of <fi(Z,Q) at 
Z = — f2, and the point Z = |f2 where <f> is equal to its 
local maximum. In the interval [— f2, ^f2], the function 
Q{Z,cj)) ~ /3exp[— <f>(Z, Q] for some constant /3. Writing 
Q{Z,n) = f3exp(-nZ 2 /2)[l~ Z 3 /3 + 0(Z 4 )}, we obtain 
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where (A)o defines the average of A with the Gaussian 
weight, Co(Z) = [3ex-p(-riZ 2 /2). 

Since A is negative for f2 3> 1, paths coalesce in 
this regime. When f2 = 0, on the other hand, A = 
{V/m 2 ) 1 / 3 ^ 12 5 / 6 r(5/6)/(2VF) > 0, which estab- 
lishes the existence of a non-coalescing phase for small 
f2 (at f2 = the system is conservative; recent studies^ of 
monodromy matrices for such systems obtain an expres- 
sion for A equivalent to ours at Q — 0). The transition 
between the two phases occurs at tt c w 0.827. The two 
examples shown in Fig. Q]arc at SI = 0.5 and SI = 1.25 
respectively. The damping which causes the maximum 
rate of coalescence is determined by the minimum of 
#i(fi)/#o(ft): this is at Q w 2.035. 

Rate of crossing of caustics. The magnitude of the flux 
determines the rate at which caustics appear on a given 
trajectory: \ J\ = (P/m 2 ) 1 / 3 /* (^), shown in Fig.03. |J| 
is largest at O = 0, j J| = {V/m 2 ) 1 / 3 T(5/6) 12 5 / 6 /(8tt 3 / 2 ), 
of the same order at Sl c , and quickly drops to zero at 
large SI. In Fig.^i caustics are still discernible. 

Statistics of the particle density. In the overdamped 
limit [v 3> 1), the momentum is approximately p(t) = 



f(x(t),t)/j, so that eq. approximated by 
dx{t) _ 1 



dt 



777,7 



-f(x(t),t) 



(14) 



In this regime, a more concrete and complete understand- 
ing of the path-coalescence effect is obtained by consid- 
ering the statistics of the density of particles g(x, t) = 
Y] j 8(xi(t) — x) . Translational invariance implies that an 
initially uniform density remains uniform, (g(x,t)) = Qq. 
Path coalescence is revealed by the density-density corre- 
lation function K,(x, x';t) = (g(x, t) g{x' , t)} — g S(x — x'). 
Because of translational invariance, the correlation func- 
tion K, is a function of Ax = x — x' only: we write 
K.(x,x';t) = A'(Ax,i). A tendency for particles to clus- 
ter is demonstrated by K(Ax, t) becoming large for Ax 
small, in the limit t — > 00. 

When x ^ I? we find that the correlation function 
satisfies a Fokkcr-Planck equation^ 



dK(Ax,t) d 2 



dt dAx 2 
The diffusion constant 
1 



[D(Ax)K{Ax,t)] 



D(Ax) 



m 2r y 2 



dt [c(0,i) - c{Ax,t)] 



(15) 
(16) 



approaches zero quadratically at the origin: D(Ax) ~ 
kAx 2 for Ax <?C £. When Ax ^ £ the diffusion constant 
approaches a constant value Dq. 

Now consider the properties of solutions of cq. 1)15(1. 
We note that eq. (|15(l is in the form of a continuity equa- 
tion, dK/dt + dj/dAx = 0, so that the integral of the 
correlation function over all Ax is a conserved quantity. 
The flux of the correlation function passing the separa- 
tion parameter Ax at time t is 

j(Ax, t) = --^-[D{Ax)K{Ax, t)] . (17) 

Consider an initially uniform distribution of density, with 
value go (corresponding to if (Ax, 0) = g^). For Ax <C £ 
the diffusion constant is an increasing function of Ax. 
Together with l|17fl this implies an initial flux of correla- 
tion towards Ax — 0. At large times, K(Ax, t) is thus 
sharply peaked at the origin. For Ax ^> £, on the other 
hand, the approximate solution of (f 1 51) is 



K{Ax,t) 



(18) 



Using the fact that K(Ax, t) satisfies a conservation law, 
we deduce that the average number of particles condens- 
ing into each trail at time t is 



N(t) ~ -=ea< 

\/7T 



(19) 



When t/r is large, and Ax <C £ (but not too close to zero) 
the flux j(A, t) is found to be approximately uniform. 
This implies 



K(Ax,t) 




7T AxVt 



(20) 




FIG. 3: Statistics of g(x, t) in the coalescing phase, a Den- 
sity correlation K{Ax,i) of the process 12H : numerical ex- 
periments (o), limiting theoretical forms 11811 . full line, and 
dashed line. Parameter values: e 2 « 1.25 x 10~ 8 , 
£ w 6.4 x 10~ 3 , and ( = 5x 10 5 r. b Numerical results verify- 
ing that K(0, t) = exp(2Kt) for the same parameter values, c 
Mean cluster size N(t), of the process (1211 : numerical experi- 
ments (o), theory 11911 . line. Parameter values: e 2 ~ 2 x 1CP 7 , 
£ » 1.6 x 10" 3 . Particles are considered to be part of a cluster 
of N if their positions are within an interval of length £. 



The 1 /A-divcrgence of (|20f) is non-integrable, so that this 
expression must fail near the origin. An exact calculation 
shows that K(0,t) = exp(2frf). 

Tracking-minima regime. When ^ > 1, there is an 
overdamped phase in which the particles follow minima 
of the random potential obtained by integrating the force 
with respect to x. A particle executes a rapid jump to a 
lower minimum whenever the minimum that it has been 
tracking disappears. 

Discrete model. The simplest model exhibiting the 
path-coalescence effect is a discrete-time random walk 
which approximates (|14ll : 



Xi(t + t) = Xi(t) + F n { Xi (t)) 



(21) 



with t = nr. By analogy with J2J we take (F n (x)) = 
and (F n (x)F n /(x')) = 5 nn <C(x — x'). Deutsche gave a 
detailed analysis of (112 1 1) . Here we remark that 1211) allows 
us to obtain the simplest possible understanding of why 
coalescence occurs. Consider the linearisation of (I2H : 

A = r 1 (log \6x(t)\/\Sx(0)\) = (log |1 + F' n \)/T . (22) 

When the magnitudes of the derivatives F' n = dF ra /dx 
are small compared to unity, Taylor expansion of the log- 
arithm gives A ~ ~\{F^ 2 ) /t, so that A is negative. This 
approximation is equivalent to 1|13|) ■ This argument indi- 
cates that the coalescence is a second-order effect. Also, 
if the random displacements are larger than their cor- 
relation length, the coalescence effect disappears: if F n 
has a Gaussian distribution, (|22f> shows that A becomes 
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positive when (F^ 2 ) exceeds 2.421 . . .. Fig. [3] shows nu- 
merical results confirming eqs. Q18|). and (|2*U|) . using 
numerical simulations of the mapping (|21|) . 

Applications. In the remainder we discuss a number of 
possible realisations of the path-coalescence effect. 

One example is the motion of liquid droplets on a sur- 
face, moving in one direction under a constant force (rain 
blown off a perspex windshield is an example of this situ- 
ation). If the surface is randomly contaminated, the wet- 
ting angle will be different on opposite sides of each drop, 
and the trajectory of the droplet will be randomly de- 
flected. We assume the surface contaminants are smeared 
over an area large compared to that of the droplets (per- 
haps resulting from cleaning the windshield with a waxy 
polish), so that nearby droplets are deflected in the same 
direction. There is no interaction between the drops un- 
less they are close enough to combine due to surface ten- 
sion: we stress that the coalescence is that of the paths 
taken by different drops, not of the droplets themselves. 
We model the motion of a droplet across the surface by a 
particle of mass m. At position r = (x, y) on the surface, 
the drop is subject to a force F(r) + Fq j, where Fq is the 
magnitude of a steady force acting in the direction of the 
unit vector j defining the y-axis and F(r) is a homoge- 
neous and isotropic random force with correlation length 
£. We assume that the particles are subjected to a vis- 
cous resistive force proportional to their velocity across 
the surface, such that the equations of motion are 

dr do 

m— =p, =F j + F(r)- 7 p (23) 

(where p is the momentum of the drop) . When the fluc- 
tuating force is weak, the trajectories are locally approx- 
imated by straight lines, with x approximately constant 
and with y increasing at a rate v y = Fq/ttij. The equa- 
tion of motion in the direction transverse to the con- 
stant force is then in the form of equation with 
f(x, t) = F x (x, v y t). In the case where the motion of the 
droplets is sufficiently damped, the trajectories coalesce 
onto fixed trails. 

The path-coalescence effect may also be relevant to 
the motion of energetic electrons in disordered solids. 
For example, the effect may be relevant to a 'branch- 
ing' observed in the flow of electrons away from a con- 
striction in a two-dimensional electron gas with very low 
scattering^. The experiment shows regions of markedly 
increased current density persisting to some distance 



from the constriction. This was explained by showing 
similarity to simulations of independent electron motion 
in the smoothly varying random potential of the doping 
atoms. This model is essentially (|23|l . with 7 = 0. The- 
oretical discussion of this systemZ has emphasised that 
caustics are important in understanding the empirical 
results. Wc remark that these experiments might show 
an even more pronounced effect if dissipation were in- 
troduced. Fig. nji is similar to the flows discussed in&, 
showing the branching effect and caustics. Adding dissi- 
pation to the electron motion would cause the paths of 
the electrons to coalesce, as in Fig. We note that 
dissipation of the electron motion could be increased by 
increasing the temperature of the system. As well as giv- 
ing a criterion for the phase transition, our expression 
for J also gives a quantitative prediction for the rate of 
formation of caustics along a given trajectory. 

There are also potential applications of the over- 
damped equations, l|14|) or 121|) . in the biological sciences, 
involving the movement of organisms in response to small 
random fluctuations in their environment. The model 
provides a mechanism through which large numbers of 
organisms can congregate without communicating. One 
example of this type is the migration of animals across a 
nearly homogeneous smooth terrain. Thus Fig. could 
be thought of as a map showing paths of animals on an 
Eastward migration. The paths of the animals will be 
deflected by small random fluctuations of topography or 
vegetation: the animals can be drawn together onto the 
same paths, even if there is no communication between 
them, and no gross features in the terrain favouring par- 
ticular routes. A second example applies to simple or- 
ganisms such as plankton which can move in response to 
changes in their environment, such as nutrient concentra- 
tion. In cases where there are small, spatially correlated 
random fluctuations of the nutrient concentration, the 
path-coalescence effect could lead to unexpectedly large 
concentrations of organisms. Such a mechanism could 
have been utilised by evolution, enabling simple organ- 
isms which cannot communicate directly to congregate 
for sexual reproduction. 
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